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(Dated:) 

Measurements are presented of the complex dynamic Young's modulus of NdNiOa and 
Ndo.esEuo.asNiOs through the Metal-Insulator Transition (MIT). On cooling, the modulus presents 
a narrow dip at the MIT followed by an abrupt stiffening of ~ 6%. The anomaly is reproducible be- 
tween cooling and heating in Ndo.esEuo.ssNiOs but only appears as a slow stiffening during cooling 
in undoped NdNiOs, conformingly with the fact that the MIT in RNiOs changes from strongly first 
order to second order when the mean R size is decreased. The elastic anomaly seems not to be associ- 
ated with the antiferromagnetic transition, which is distinct from the MIT in Ndo.esEuo.ssNiOs. It is 
concluded that the steplike stiffening is due to the disappearance or freezing of dynamic Jahn- Teller 
(JT) distortions through the MIT, where the JT active Ni'^^ is disproportionated into alternating 
Ni^^* and Ni"^~*. The fluctuating octahedral JT distortion necessary to justify the observed jump 
in the elastic modulus is estimated as ~ 3%, but does not have a role in determining the MIT, since 
the otherwise expected precursor softening is not observed. 



I. INTRODUCTION 

There is still uncertainty on the exact nature of 
the metal-insulator transition (MIT) in the i?Ni03 
perovskitesji"— Following a systematic study with sev- 
eral R ion sizeSfi the MIT in these perovskites has been 
associated with the opening of a gap between the O 2p 
band and the Ni 3d upper Hubbard band when the R 
ion size and/or temperature are decreasedi^i^ In fact, the 
i?— O bond is too short with respect to the Ni— O one 
for the cubic perovskite structure and the mismatch is 
relieved by tilting of the NiOg octahedra. Both a smaller 
R size and cooling enhance the mismatch and hence tilt- 
ing, for steric reasons and because of the larger thermal 
expansion of the R—0 bond. This results in a further 
reduction of the angle between Ni— O— Ni bonds with re- 
spect to 180°, reduces the overlap and hence the width of 
the O 2p and Ni 3c? bands and finally opens a gap between 
them, causing the MIT.— Indeed, with increasing the R 
ion size the temperature Tim at which the MIT occurs 
decreases, or the metallic phase becomes more stable.^ 
In this manner it is possible to rationalize the phase di- 
agram of i?,Ni03 of T vs the R ion size, and a similar 
behavior is found in cobaltites.^ 

On the other hand, Ni'^^ has an electronic configura- 
tion similar to that of Mn'^+ , with filled t2g triplet states 
(thanks to Hund's rule in the case of Mn) and one elec- 
tron in the eg doublet, whose degeneracy can be lifted by 
tetragonal and orthorhombic distortions according to the 
Jahn- Teller (JT) effect. The physics of Mn perovskites 
is dominated by the Jahn- Teller coupling between these 
eg electronic states and the distortions of the octahe- 
dra, which is a cause of orbital ordering (00), and the 
same might be expected for nickelates. On the contrary, 
according to the early diffraction experiments the octa- 



hedral distortions in nickelates are extremely small or 
nuU.^— Later, starting with the nickelates with smaller R, 
it has been found that the MIT is accompanied by a sub- 
tic orthorhombic to monoclinic structural change, with 
charge disproportionation (CD) 2m^+ Ni^^^ 
and charge ordering (CO) into alternately expanded and 
contracted NiOe octahedra along the three directions?^ 
and some octahedral distortion of the larger Ni^~'^06 oc- 
tahedra. A similar result has been recently found also 
with the larger R ions Nd^ and Dy.l° Such JT distortions 
are found only in the larger and more ionic Ni'^^'^Og oc- 
tahedra and are at least one order of magnitude smaller 
than those observed in manganites. 

Notice that in a naive picture CO and 00 should not 
occur in the same set of Ni ions, since Ni'^"'" is JT ac- 
tive, while Ni^+ and Ni^"*" are not, but to what extent 
CO and 00 occur in nickelates is still subject of con- 
troversy. Just because of the subtleness of the structural 
changes attributable to CO or 00, recourse has been had 
to resonant X-ray scattering. --^''4, These experiments in 
NdNiOs are interpreted as evidence of partial CO rather 
than 00 occurring at the MIT. According to this view, 
the degeneracy of the eg orbitals would be lifted already 
in the metallic state, where the 00 would be reflected by 
the tilts of the octahedra; the MIT would be due only to 
CD, accompanied by a small change of tilt angle It is 
debated to what extent these experiments may provide 
information on the contribution of CO and 00 to MIT^i^ 
but the view that MIT is due to CO rather than 00 
received further support from theory and experiments 
under hydrostatic pressure in LuNiOs,^^ and from the 
above mentioned diffraction studies on R — Nd and Eu, 
where alternated small and large octahedra are found.— 
Nonetheless, recent experiments under pressureii suggest 
that the previous ones^ were affected by an unwanted 
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uniaxial component of pressure, and exclude a role of CO 
in the MITj^ rather, dynamic JT deformations would be 
responsible for large isotope effects in Tim-^^ 

Additional experiments that have been performed in 
order to probe CD with CO and 00 at the MIT in 
nickelates include Mossbauer— and x-ray absorption 
spectroscopies.— Again, there is no clear picture, with 
possibility of 00 ^^'^° or C022r— even in the metallic 
stated 

A similarly confused situation exists for NaNi02 and 
LiNi02 , which are also composed of NiOe octahedra but 
sharing the edges instead of corners While NaNi02 
undergoes a phase transformation with clear coopera- 
tive JT distortion, LiNi02 has little or even reverse JT 
distortion,-^ which is not yet explained. Some of the pro- 
posed explanations involve the peculiar geometry of the 
LiNi02 lattice, where the triangular Ni sublattice frus- 
trates AFM correlations, but other possible causes may 
apply also to iJNiOa, like the fact that the JT distortion 
exists but is incoherent2£ or with very short coherence 
lengthf2^ or is hindered by the electron delocalisationJ^ 

From the above summary and recent reviews^ it ap- 
pears that, in spite of the wealth of experimental data 
on the i?Ni03 series, the presence and role of CD and 
CO and/or 00 at the MIT is not clear. In the Mn 
perovskites important and quantitative information has 
been obtained on such issues by studying the associated 
elastic anomaliesr^i^ in Ni perovskites no such type of 
investigation has ever been reported, except for a prelim- 
inary report of the present results. — Here anelastic spec- 
troscopy measurements on Ndi-ajEu^jNiOa are presented 
and discussed, which shed new light on these issues. 



II. EXPERIMENTAL 

Polycrystalline samples of Ndi-ajEu^^NiOs {x — Q and 
0.35) were prepared from sol-gel precursors, sintered at 
T = 1000 °C, and under oxygen pressure of 80 bar. De- 
tails of the synthesis process for preparing these samples 
are described in details elsewhere. All samples were 
characterized by X-ray powder diffraction in a Brucker 
D8 Advance diffractometer. The x-ray diffraction pat- 
terns showed no extra reflections due to impurity phases, 
and indicated that all samples have a high degree of crys- 
tallinity. The samples were two bars with a; = (NENO) 
and with x = 0.35 (NEN35), both with approximate size 
43 X 6 X 0.7 mm'^ and a density p ~ 2.9(3.1) g/cnr^ for 
NEN0(35), about 40% of the theoretical density. Such a 
high porosity makes the evaluation of the absolute value 
of the complex Young's modulus E = E' + iE" problem- 
atic, but did not affect the quality of the measurement of 
its temperature dependence. This was measured by elec- 
trostatically exciting the flexural vibrations of the bars, 
which were suspended on thin thermocouple wires in vac- 
uum, and with one face made well conducting also below 
the MIT with Ag paint. It was possible to excite various 
flexural modes, whose resonance frequencies are ideally 



in the ratios /if : /2F : /3F : /sf = 1 : 2.800 : 5.487 : 
13.55. The fundamental frequency is given by^ 

Af = 1.028;i^ (1) 

where / is the length, h the thickness, p the density of 
the bar; at room temperature it was /if 0.6 — 1 kHz. 

In the Discussion we will need an estimate of the abso- 
lute value of the Young's modulus, which is found with 
Eq. ([T]). The effective Young's moduli at room temper- 
ature, uncorrected for porosity, were E = 4.5 GPa for 
NENO and 12.5 GPa for NEN35, which are 5-10 times 
smaller than the typical values E ^ 100 GPa found in 
ceramics of similar type. Sources of error are the sam- 
ple porosity, deviations from the shape of a rectangular 
parallelepiped, low homogeneity and the presence of the 
Ag conductive layer, wrong identification of the vibration 
modes and mixing of flexural with torsional modes. The 
largest source of error in our case seems to be the ex- 
tremely high porosity, p ~ 0.60. The other causes should 
be of minor importance; for example, the NENO sample 
had the lowest value of E, but its ratios /sf/Af = 5.60 
and /5F//1F — 13.66 were in excellent agreement with 
the theoretical values given above. The correction for 
such high porosities are very unreliable; to give an idea, 
the Young's modulus of porous ceramics are corrected 
with empirical expressions like (1 — p)" with n ~ 1,— or 
exp {—bp) with 2.1 < 5 < 2.8.'^'' These expressions would 
enhance E in our case by a factor between 2.5 and 5.2. 
While this brings E of NEN35 to reasonable values up to 
66 GPa, the value of NENO would remain below 25 GPa. 
The large difference between the two may in part be ac- 
counted for by different types of porosities, in terms of 
shape and connectivity of the pores, and in part due to 
the Eu substitution. In the discussion we will assume 
S - 65 GPa. 

The anelastic spectra are displayed in terms of the real 
part of the Young's modulus or its reciprocal, the com- 
pliance s = E^^ = s' — is" , and of the elastic energy loss 
coefficient Q^^ = s" / s' . 

III. RESULTS 

In Fig. [T] is reported the anelastic spectrum of NENO 
measured during cooling (empty symbols) and subse- 
quent heating (filled symbols) measured exciting the first 
odd flexural modes at 0.63, 3.5 and 8.6 kHz during the 
same run. During cooling, below 190 K there is a pro- 
gressive stiffening of the Young's modulus in excess of 
the slight linear anharmonic stiffening observed at higher 
temperature; such a stiffening proceeds until the lowest 
temperature reached by us, although in the end it starts 
levelling off. On heating, the modulus resumes the al- 
most linear and weak temperature dependence except for 
a steplike softening followed by a dip at 191 K. It is clear 
by comparison with resistivity, specific heat and magne- 
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FIG. 1: Young's modulus E and elastic energy loss coefficient 
of NENO measured on cooling (empty symbols) and heat- 
ing (filled symbols) at 0.63 kHz (black), 3.5 kHz (green) and 
8.6 kHz (red) (colors online). 



tization measurements that the steplike change of E (T) 
with large temperature hysteresis is due to the MITj^i^i 

The (T) is perfectly reproducible on heating and 
cooling, except for the sharp peak in correspondence with 
the cusped softening, both of which appear only during 
heating. In addition, there are other three peaks: PI at 
230 - 240 K, P2 at 171 K and P3 at 100 - 110 K. While 
the temperature of P2 is independent of frequency, the 
other two increase with frequency, and therefore indicate 
thermally activated processes whose maxima occur when 
the condition ujt ~ 1 is met,— where uijl-K is frequency 
and r (T) a relaxation time. We will ignore these repro- 
ducible but small peaks. 

The spectrum of NEN35 (Fig. [2]) is similar to that 
of NENO during heating, but the MIT occurs at Tim = 
288 K, if identified with the cusp in the real part, and 
the hysteresis is completely absent, the only difference 
between heating and cooling being the intensity of the 
peak in real and imaginary parts at the MIT. The curves 
measured at a frequency five times higher are not shown 
for clarity, since they are identical except that the peak 
in Q"^ is ^ 15% higher, whereas in NENO it is nearly 
identical. 

The other small peaks in the are absent in NEN35. 
In Fig. [2] are also reported the curves measured on 
NdNiOs during heating, rescaled in temperature by a 
factor 1.504. The rescaled elastic anomaly of NdNiOa 
is slightly sharper and narrower in the real part, with a 
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FIG. 2: Young's modulus E and elastic energy loss coeffi- 
cient of NEN35 measured on cooling (empty symbols) 
and heating (filled symbols). For comparison, the anelastic 
spectrum of NENO on heating is reported (line) after scaling 
its tempearure to let its Tim coincide with that of NEN35 
(colors online). 



smaller dip, and almost coincident in the absorption. In 
NEN35 the antiferromagnetic (AFM) transition occurs 
at Tn = 240 K < Tim,?^ but there is no clear sign of it 
in the anelastic spectra. 



IV. DISCUSSION 

The overall phenomenology of the anelastic anoma- 
lies found in Ndi^ajEuajNiOa appears to be mainly re- 
lated to the MIT and in agreement with what is known 
from other techniques. In NdNiOa the MIT occurs at 
Tim = 7n — 192 K and has marked first order charac- 
ter, with a large hysteresis between cooling and heating. 
When decreasing the mean i? ion size by Eu doping, the 
MIT shifts to higher temperature and becomes second 
order, while Tn shifts below Tim. The anelastic anoma- 
lies also appear to be marginally affected by the mag- 
netic transition, since the heating curves are very similar 
to each other regardless of the fact that Tn coincides 
with Tim or not. The present data can be compared 
with those of resistivity, magnetization and specific heat 
taken on the same materials . ^^'"^^ The hysteresis is due 
to the fraction of material that remains metallic in the 
insulating phase, and the fact that the step in p (T) (see 
Fig. 7 in Ref. 1341 1 is narrower than that in E (T) during 
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cooling is explained by the fact that the latter probes 
the bulk fraction of insulating and metallic phases, while 
resistivity rather probes percolating conductive paths. 



A. Comparison with the specific heat anomaly 

It seems that the dip in the modulus at Tim and 
the stiffening below it have different origins, since with 
X = 0.35 the dip and accompanying peak in differ in 
amplitude between heating and cooling, while the step in 
E is reproducible. This observation allows us to separate 
the two anomalies and compare the dip in the modulus, 
or peak in compliance s = , with the peak in specific 
heat. The curves 1 and 2 in Fig. [3] are the normalized 
compliance on heating and cooling, after subtracting a 
linear background Sbg] the reference so is the extrapola- 
tion of Sbg to K. Besides the peak at 290 K and the 
step below that temperature, there is an additional shal- 
low rise of the compliance below 250 K. The shape of 
the latter is very uncertain, due to the background sub- 
traction. The anharmonic phonon contribution to the 
temperature dependence of the elastic constants is gen- 
erally almost linear and levels off at low temperature, but 
being unable to distinguish the anharmonic contribution 
from the other anomalies in the temperature dependence 
of s, we chose to use a simple linear interpolation over the 
whole 60 — 350 K range. Curve 3 is the difference between 
curves 1 and 2 and it has been subtracted from them af- 
ter multiplication by scale factors in order to remove the 
peaked component. It turns out that the best values of 
the factors are 2.00 and 1.00. The resulting curves 4 and 
5 are coincident and have only the step below 290 K plus 
the shallow rise. 

Both the peak Spcak (curve 3) and the excess compli- 
ance below 250 K can be put in relationship with the ex- 
cess specific heat Cex measured on a similar sampled* In 
fact, a softening proportional to the specific heat anomaly 
is expected at a magnetic transition^SS, and also at other 
types of transitions under the rather general condition 
that their critical temperature depends on pressure or 
more generally on stress.^ A possible identification of the 
two anomalies in Cex is the peak at ~ 300 K with the CO 
and/or orbital transition concomitant with the MIT, pos- 
sibly with participation of spin degrees of freedom, '^^ and 
the broad contribution below 250 K with the onset of long 
range AFM order. The latter is expected to produce a 
jump AS" = i?ln 2 in the entropy of the S = ^ spinsr^i^S 
but the jumps of the entropy S = /dT Cex/T are only 
AS" = 0.15i? between 130 and 260 K and AS* ~ 0.12i? 
between 260 and 320 K. The fact that the actual step is 
smaller than expected from complete AFM ordering has 
already been noted in other iJNiOa nickelates, and has 
been interpreted as due to AFM correlations also above 
Tn, so that only a relatively small loss of entropy is in- 
volved in the long range AFM ordering i^^'^^ 
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FIG. 3: (a) Curves 1 and 2: compliance s = E'^ of NEN35 af- 
ter subtraction of the anharmonic background Sbg, measured 
on heating and cooling. Curve 3 is the difference between 
them, Speak- Curve 4,5 are curves 1,2 after subtraction of the 
peaked component 3. (b) Excess specific heat extracted from 
Ref . [2i. 



B. Step-like stiffening below the MIT 

The main feature of the elastic anomaly at both com- 
positions is the stiffening below the MIT, whose step- 
like nature is made evident after the subtraction of the 
peaked component in curves 4 and 5 of Fig. [3) This 
is quite an unusual observation, since phase transitions 
are generally accompanied by softening rather than stiff- 
ening. As discussed below, if some strain component is 
coupled linearly with the order parameter of the transi- 
tion, then its elastic constant presents a cusped softening 
at the critical temperature and is therefore followed by 
restiffening to the background elastic constant on further 
cooling. This is not the case of curves 4 and 5, which lack 
any sign of precursor softening divergent near Tim ■ 

There are several elastic studies of the CO and 00 
transitions in perovskite manganites and in other com- 
pounds, but most of them concentrate on the precursor 
softening above the transitions. In few cases the sub- 
sequent restiffening below the transitions has been con- 
sidered, especially if it appears of excessive sharpness 
and amplitude with respect to the precursor softening. 
In these cases recourse has been had to different sets 
of fitting parameters above and below Tco ,— or to a 
Landau free energy expansion with coefficients unrelated 
to CO or 004i. In Ndo.sSro.sMnOs a steplike stiffen- 
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ing at the MIT has been interpreted as the effect of the 
renormahsation of the elastic constant by the conduc- 

tion electrons with a factor (l + g^x [T)) , where g is 
the electron-phonon coupling constant and x C^) is the 
electron susceptibility, identified with the magnetic sus- 
ceptibility Xm'^ In that case, the MIT coincides with a 
transition from FM to AFM and Xm has a negative step, 
which becomes positive in the renormahsation factor of 
the elastic constant. This type of analysis is not appro- 
priate to Ndi_j^Eua;Ni03, whose magnetic susceptibility 
has an almost imperceptible decrease below Tn, super- 
imposed to a stronger rise with cooling, even after sub- 
traction of the contribution of the Nd and Eu ions i^'*'^'^ 
Lacking an adequate precedent for a satisfactory inter- 
pretation of the step component of the elastic compliance 
of Ndi-ajEu^^NiOa, we first review what kind of elastic re- 
sponse is expected at a MIT or magnetic transition, and 
then propose an interpretation based on the response of 
the JT distortions, from a point of view different from 
that usually adopted. Detailed discussions of the elastic 
anomalies expected at various types of phase transitions 
can be found in many review articles and books, usually 
based on the Landau theory with the coupling between 
strain and order parameter included;^"— and therefore 
we only quote what is strictly necessary to our discussion. 



C. Landau analysis 



D. Magnetic transition 

Magnetic transitions cause elastic anomalies through 
essentially two mechanisms: exchange striction and mag- 
netostriction. The first mechanism causes cusp-like soft- 
ening in the elastic constants involving the strains that 
change the atomic distances and hence the exchange con- 
stant. The magnitude of such anomalies is of the order 
of Even smaller step-like softening can be caused 

by magnetostrictive coupling with the lattice, namely 
by a term —BsijSiSj linear in strain and quadratic 
in spin variable or magnetization, which, according to 
Landau's theory,^^ produces below the magnetic tran- 
sition a negative step in the modulus M of magnitude 
—2B^/M. Such a softening is usually observed at mag- 
netic transitions j^^i^i^ 

On the other hand, in the FM transition of YTiOa, a 
stiffening is found, whose origin has not been explained, 
but also involves 00. In the hexagonal quasi- ID AFM 
multiferroic YMnOa, there is 1% stiffening below Tn in 
Cii and 3% in Cee due to biquadratic coupling^ (see 
previous paragraph). In fact, below Tn = 75 K the tem- 
perature dependence of ACu (T) / Cu follows the squared 
OP, 5*2 (X (1 - T/T^f^, downto 2 K, which appears dif- 
ferent from the abrupt rise in Ndi-^^Euj^NiOs. In the 
present case, the involvement of magnetic effects in the 
steplike stiffening is made even more unlikely by the fact 
that in Ndo.esEuo.ssNiOs the magnetic transition occurs 
at Tn < Tim- 



Near a phase transition the free energy can be ex- 
panded in powers of the order parameter(s) Q. In the 
present case the order parameter (OP) can be one or 
more of the symmetrized charge fluctuation coordinates 
for describing charge ordering, or the quadrupolar or- 
bital operators of the Ni ions for orbital ordering, or the 
staggered magnetization for the AFM transition. By in- 
cluding powers of Q up to the 6th order, it is possible 
to reproduce both first and second order transitions, and 
including coupling terms containing both Q and strain e, 
it is possible to deduce the effect of the relaxation of the 
OP under stress on the elastic constants. In the case that 
the OP is strain itself or is linearly coupled to it with a 
term \eQ, the elastic constant coupled to the OP ideally 
vanishes at the transition temperature Tq , or at least has 
a negative cusp. If bilinear coupling is forbidden because 
no strain has the same symmetry of the OP, the next 
coupling term, ^eQ^, causes a negative step at Tc in the 
the elastic constant coupled to the OP, possibly includ- 
ing a weak restiffening below Tq. A biquadratic coupling 
ue^Q^ adds to the elastic constant a contribution 2v (Q^) 
below Tc, which can be a linear or saturating rise or de- 
crease, depending on the temperature dependence of the 
OP and on the sign of the coupling constant i/. Of all 
these terms, the latter is the only one that might pro- 
duce a stiffening below Tc without precursor softening; 
other terms are possible but less important in the great 
majority of cases. 



E. Cooperative Jahn- Teller phase transition and 
orbital ordering 

Following the notation of Hazama and coworkers 
the quadrupolar moments of the Cg orbitals of Ni'^"'" can 
couple with tetragonal e" = {2ezz — £xx — £yy) /V^ and 
orthorhombic = Exx — £yy strains (both with the same 
Eg elastic constant C" = (Cu — C12) /2), and interact 
with each other causing the cooperative JT transition:^ 
The relevant part of the Hamiltonian of N Ni'^+ ions per 
unit volume, referred to the unit cell volume uq, is 

H = -voN 9'fO-<e^ " ^oN ^ g' (O^) O^,. (2) 

'Y—u^v j—u^v 

where the quadrupolar operators O2 = 
[2II — — ly) /-y/S and O2 — 11 — ly correspond, apart 

from a geometrical factor a — {(j)eg \ Oeg \4'eg) = 2-\/3, to 
the occupation numbers of octahedra with tetragonal 
and orthorhombic JT distortions refiecting the symmetry 
of their Sd orbitals. The first term is the linear coupling 
to the external strain, with a coupling strength and 
causes softening of the C" elastic constant, while the 
second term is the elastic interaction among the orbitals 
in the mean field approximation with strength g' , and 
determines the type of 00. The resulting softening of 
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the C elastic constant is 



1 + 9'X{T) 



(3) 



where Cq is the background elastic constant and the sus- 
ceptibility 



xiT) 



(4) 



is proportional to the contribution to the compliance s = 
from noninteracting orbitals, As (T) = (^) x (T)- 
Equation ([3]) can be rewritten as 



T- 



Too \ 



(5) 



where Q = a g'/Zce , Too = 6 + Tjt and Tjt = 
N a? / [kBCo) . In this form it is clear that the elas- 
tic constant vanishes at Too, which is the onset of the 
00 or cooperative JT transition, and that stiffens again 
on further cooling. At higher temperatures, when T ^ Q 
and hence g'x (T) 1, Eq. (jS]) shows that the softening 
induced by the orbitals freely responding to the applied 
stress is simply given by Eq. Q and is cx 1/T. 



F. Charge ordering 

The charge fluctuations can be expanded into fluctu- 
ation modes Qj acting as OP of the CO transition, one 
or some of which may be linearly coupled to the strain 
e-y having the same symmetry with terms —g^QjEj. As 
also shown in Appendix B, if such terms are included in 
the expansion of the free energy in powers of up to 
the fourth order, a second order transition at a temper- 
ature Tco results, and the renormalized elastic constant 
can be written exactly as in the case of the 00 tran- 
sition, Eq. ([5]), where O is the temperature at which 
the term ^ (T — 8) of the Landau expansion van- 
ishes and Tco = O + AT, AT = g^/ (aoC^), with 
the elastic constant appropriate to e^, is the tempera- 
ture at which the CO transition occurs with onset of 
spontaneous strain of type e-y. Therefore, 00 between 
two orbitals and CO from condensation of a single charge 
fluctuation mode produce the same precursor softening 
on the appropriate elastic constant (= C" for 00 
within the eg doublet). The restiffening below the tran- 
sition is discussed in Appendix B. Examples of materials 
similar to iJNiOa where both transitions have been recog- 
nized in ultrasonic experiments are Lai_a;Sra;Mn03 with 
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oo 



and Pri_j;Ca2,Mn03 with Tco > Tqo^ 



active orbitals or between ionic charges would appear in 
the coupled elastic constant as a clear precursor soften- 
ing that diverges at the transition, followed by restiffen- 
ing on further cooling. This is clearly not the case of 
Ndi_a;Eu2,Ni03, since the Young's modulus presents a 
weak and almost linear stiffening downto a temperature 
very close to the MIT, which excludes any mechanism 
coupled to strain. This is also shown quantitatively in 
Sect. I. 

The picture we propose is the following. The MIT is an 
electronic transition, presumably occurring at the open- 
ing of the charge transfer gap during cooling)^ when 
the R—0 bonds, having larger thermal expansion of the 
shortest Ni— O bonds, shorten faster than the latter and 
enhance the buckling of the Ni— O— Ni bonds beyond a 
critical threshold. The driving mechanism for this is the 
thermal contraction due to lattice anharmonicity, and 
there is no effect on the elastic constants other than the 
usual anharmonic stiffening. In the metallic phase the 
octahedra would be JT distorted with a fluctuation rate 
much faster than our measuring frequency, so that their 
contribution to the lattice softening would be noncritical, 
namely Eq. ^ with g' = 0, and the 00 transition would 
occur well below Tim. The MIT would inhibit such a 
softening, possibly due to the concomitant CO that trans- 
forms JT active Ni'^'^ into (partially) inactive Ni^^^ , or 
making the distortions static, or both. This would make 
the insulating phase stiffer than the metallic one, so that 
the disappearance of the metallic phase would be accom- 
panied by elastic stiffening. 

In what follows we will evaluate the JT distortion nec- 
essary to explain the steplike stiffening in Figs. [1] and 
[2] with the above mechanism and test the consistency of 
this explanation with the available data. In order to re- 
late the amount of elastic softening to the JT distortion 
of the octahedra, rather than to the coupling coefficient 
g, it is convenient to recast the above results in the usual 
formalism of anelastic relaxation from point defects.— 
It may seem unrealistic to describe distorted octahedra 
sharing the O atoms, and hence strongly interacting, as 
almost independent defects carrying a distortion, but this 
is equivalent to the formulation of Sect. IIVE[ namely the 
00 approach to the cooperative JT effect^S. In the case 
of nickelates it is certainly justified by the smallness of 
the JT distortions, which make uncorrelated fluctuations 
less energetically unfavorable. If we neglect the interac- 
tion term cx g' , the hamiltonian ([2]) is equivalent to that 
of a system of independent defects of type 7 = u, w with 
double-force tensor p'^ , whose interaction energy with an 
external strain e is^ 



7 u 



G. Softening in terms of octahedral JT distortion 
and its disappearance below the MIT 

Prom the above discussion it appears that a transi- 
tion whose driving force is the interaction between JT 



where with c-y — vqN are the molar concentra- 
tions of defects, and the summation over the cartesian 
components i,j selects the strain component of the same 
symmetry of the defect, . It appears therefore that 
NgjOy in Eq. ^ is equivalent to Cjp"^ , although the 
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separate correspondences NO^ — >■ c-y and — >■ p'' do not 
necessarily hold. In turn, the double force tensor can be 
expressed in terms of the elastic quadrupole Aj^ (usually 
called elastic dipole-^) as pj^ =J2ki ^ijki^kn where C is 
the fourth-rank elastic stiffness tensor. The tensor is 
adimensional and is the strain of a unity concentration of 
defects of type 7, in the present case the tetragonal or or- 
thorhombic strain of the octahedra, so that the anelastic 
strain due to defects is 



where the c-y obey Boltzmann's statistics and depend on 
the applied stress a through the defect energies E^y — 
—vop'^e^ — ~'Vo\''a~f. For only two types of defects with 
the same symmetry and coupled to strain through the 
same symmetrized elastic constant C" the relationship 
between A and p is simply A = p/C , and the contribution 
to the compliance, s' = 1/C", ig^^ 



, , de^" (AA)' 
OS ~ — — oc c- 



da 



(6) 



where A A = |A„ — Ai,| and the numerical factor differs 
from those usually reported for reorientation of defects 
with the same symmetry, since we are dealing with the re- 
laxation between a tetragonal and an orthorhombic state 
rather the reorientation of the principal axes of A among 
equivalent crystallographic directions. 

We can carry the analogy with point defects further 
into the dynamic response. If r is the relaxation time for 
an electron to change orbital, the softening Eq. © ac- 
quires a frequency dependent factor (1 -I- icor) ^ notice- 
able when temperature is lowered enough to slow and 
make it comparable to the measuring angular frequency 
this has also been discussed in connection with JT ef- 
fect in UO2.— The concentration c of eg electrons would 
be 1 for i?Ni03 if there were no disproportionation of 
Ni'^+. The MIT may reduce the magnitude of Ss' both 
reducing c through the concomitant disproportionation 
of Ni^+ and/or making null the frequency factor if the 
distortions become quasistatic with ljt 3> 1. 

In Appendix A we define the octahedral distortions in 
terms of the symmetrized tetragonal and orthorhombic 
strains T" and with magnitude A, and evaluate the re- 
sulting polycrystalline average of the relaxation strength 
of the Young's modulus, A = SE^^ / E^^ , as 



(A) 



■E 



(7) 



H. Estimate of the Jahn- Teller distortion necessary 
to explain the stiffening below the MIT 

As a first check that this mechanism can explain the 
stiffening at the MIT, let us suppose that on passing 



across the MIT the JT distortion totally disappears to- 
gether with the Ni'^+ ions; therefore we equate the above 
(A) with c = 1 to the overall jump of 6s/ sq ~ SE/Eq. 
The largest source of error in evaluating A is the abso- 
lute value of the Young's modulus, which we assume to be 
£: - 65 GPa (see Sect. HT]), and together with A = 0.056 
for NEN35 yields SE'^ ~ IQ-^^ cmVerg; the cell vol- 
ume is^ z;o = 55 X 10-24 cm3 and Eq. ^ with T = Tim 
gives A ~ 0.031, namely 3% octahedral and tetragonal 
distortions. 

In order for this interpretation to be self consistent, it 
should turn out that the temperature Too = Q + Tjt for 
00 is smaller than Tim. There is no constraint to the 
smallness of the intersite contribution O, which can also 
be negative for antiferroquadrupolar coupling. Instead, 
the JT energy is determined by the coupling with uni- 
form strain, A, with Tjt — Na?g^ jk-QCl^, where g is de- 
termined by the equivalence Nga = cp = cC'A. The Cq 
elastic constant is unknown and we set it equal to E, find- 
ing g/fcs ^ 2300 K and Tjt = 240 K, therefore smaller 
than Tim. The same estimate for NdNiOa with A = 0.062 
and T = 190 K yields A = 0.026, g/k^ ~ 1900 K and 
Tjt = 177 K, again smaller than Tim. 

Up to now it is shown that a CD with complete sup- 
pression of the JT effect can explain the stiffening below 
the MIT, but it is also necessary to determine the ef- 
fect of an incomplete charge disproportionation 2Ni'^+ — >■ 
Ni^+^-l- Ni^-'' with 5 <1. The above treatment suggests 
a naive analogy of JT active octahedra with point de- 
fects, whose concentration is c = 1 in the conductive 
phase and becomes c = 1 — 5 in the insulating one, with 
unchanged distortion A. However, this is not the pic- 
ture proposed after the diffraction experiments, where 5 
is deduced from the change of valence of Ni'^+ , and this 
is in turn extracted from the average Ni— O bond lengths 
together with the bond valence model.— In this picture, 
all the octahedra undergo disproportionation, and are 
alternated along the three principal directions as small 
Ni'^"'"'' and large Ni'^"'' octahedra, only the latter sup- 
porting some residual JT distortioni^iiS From this point 
of view, a more sensible hypothesis is that, when enter- 
ing the insulating phase, in Eq. ([T]) c passes from 1 to 
1/2 and the JT distortion A is reduced by an amount 
depending on 5. If we suppose a linear relationship be- 
tween A and (5, then we should set A — ^ A (1 — 5) into 
Eq. ([7|), and repeating the above analysis with the es- 
timate 5 ~ 0.28 for NdNiOs^iii and c = 1/2 we would 
find A = 0.036 and Tjt = 330 K. The latter value, being 
slightly larger than Tim, seems at first inconsistent with 
the present analysis, but there are several arguments that 
this is not necessarily the case. A first point is that we 
ignored the correlations between octahedra, which in the 
perovskite structure with octahedra sharing the vertices 
is expected to be antiferroelastio^ and corresponds to a 
negative 8 compensating Tjt, as in Pri_a;Caa;Mn03.— 
Another point is that the JT distortion may not fall lin- 
early with 5 but faster. Finally, if the JT distortions 
remained frozen below the MIT, then the kinetic factor 
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1/ 1 + (wr)^ would pass from 1 to on crossing the 

MIT and the initial safe estimate of A of Eq.Q with 
c = 1 would again be correct. 

Once established that the stiffening below the MIT 
can be consistently explained by vanishing or blocking 
of dynamic JT distortions A, it is necessary to check 
the compatibility of such disordered distortions above 
Tim with the structural studies on NdNiOa. We refer 
to the most recent experiment with synchrotron powder- 
diffraction^ where the Ni— O bond lengths are reported 
to be / ~ 1.9 A, and the B-factors of the apical and equa- 
torial O atoms are 0.72 and 1.13 at room temperature, 
resulting in thermal plus disordered displacements of 0.1 
and 0.12 A respectively. The JT distortions estimated 
here of A ^ 0.03 require changes of the bond lengths of 
Sl = Xxl^ 0.06 A, which, being dynamically disordered, 
would appear only in the Debye- Waller factors, and are 
about half of the B-factor displacements at room temper- 
ature in NdNiOs. The JT distortions required to explain 
the stiffening at the MIT are therefore compatible with 
the diffraction experiments: above the MIT they are dis- 
ordered and account for half of the thermal/disordered 
displacements of the O atoms; below the MIT they are to- 
tally or partially absent due to the disproportionation of 
the JT Ni^^^ into Ni'^^^, or simply frozen. This may also 
explain EXAFS experiments indicating that the splitting 
of the Ni— O bond lengths persist also above the MIT in 
i?Ni03 with R = Pr, Nd, Eu.^^ This fact has been ex- 
plained in terms of localized CO phase also in the metallic 
state, but might be due to the disordered JT distortions. 
Another case where the anomalously large mean-square 
displacements suggest JT distortions is a recent neutron 
diffraction study on PrNiOsj^ 



I. The peaked softening at Tim 

Up to now we dealt with the steplike stiffening be- 
low the MIT, which however is accompanied by a peaked 
softening and absorption. The narrow shape of such an 
anomaly suggests that it is due to fluctuations associated 
with the MIT, but the absorption is almost independent 
of frequency, while the critical absorption is expected to 
be proportional to it 42. 

The real part of the peak, corresponding to curve 3 
in Fig. [31 is plotted in Fig. [Has Spoak + sq, cleared of 
the step and temperature dependent part of the back- 
ground. We tried to fit with the expressions appropriate 
for CO and 00, namely s = 1/ (Ccq/oo + C'bg) where 
Cco/OO {T) is given by Eq. ^ above Tco/Tqo and in 
Appendix B in the low temperature phase. The param- 
eters were Too = Too = 289 K, AT = Tjt = 14.4 K 
and Co ~ 0.02, Cbg = Sq ^ = 0-98 for both cases. The 
resulting curves, very close to each other, fit very poorly 
the experimental data. In particular, the observed peak 
is steeper in the high temperature side, while the fitting 
curves are steeper in the low temperature phase, as ex- 
pected. This discrepancy cannot be accounted for by 
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FIG. 4: Peaked component of s' (curve 3 of Fig. [3| fitted 
assuming CO or OO at 289 K. 



relaxation of domain walls in the low temperature phase, 
since relaxations depend on the measuring frequency and 
the anomaly we observe does not. The failure of repro- 
ducing the peaked component as softening from 00 or 
CO should be considered as an indication that the ori- 
gin of the peak is different and we leave it as an open 
question. 



V. CONCLUSIONS 

We measured the complex dynamic Young's modulus 
of NdNiOa and Ndo.esEuo.ssNiOa, in order to study the 
elastic anomalies associated with the MIT and magnetic 
transition. The main features are a sharp stiffening of 
~ 6% below the MIT, perfectly reproducible during cool- 
ing and heating in Ndo.esEuo.asNiOs and with a broad 
hysteresis in NdNiOs, accompanied by a narrow dip at 
the MIT. The change of the MIT between first- and 
second-order character at the two compositions agrees 
with the known phase diagram of i?Ni03 as a function 
of the R size. The differences with respect to elasticity 
measurements on Mn perovskites are the absence of pre- 
cursor softening associated with the MIT, the sharpness 
of the transition during heating or when it is second or- 
der, and the broadness of the hysteresis when it is first 
order. Both the narrow softening at the MIT and the 
stiffening below it seem not to be associated with spin 
ordering, which in the Eu-dopcd sample occurs at a tem- 
perature lower than the MIT. 

By comparing cooling and heating runs, it was possible 
to separate the steplike stiffening from the other anoma- 
lies, which present similarities with the excess specific 
heat, and whose relationship with spin and charge order 
have been briefly considered. A fit of the dip at the MIT 
in terms of CO has been made, but its quality is too poor 
to draw a conclusion. 
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In order to explain the rather unusual steplike stiffen- 
ing, possible mechanisms producing elastic anomalies at 
magnetic, charge order and orbital order transitions have 
been considered, and it is concluded that the stiffening is 
due to the disappearance or to freezing in the insulating 
phase of dynamic JT distortions, due to the charge dis- 
proportionation and ordering that accompanies the MIT 
and transforms JT active Ni'^^ into Ni'^^'^. The driving 
force for the MIT is not orbital ordering or any mecha- 
nism linearly coupled to strain, in view of the absence of 
precursor softening. This is in agreement with the idea 
that the MIT is caused by the opening of a gap when the 
orbital overlap drops below a threshold during cooling. 
In fact, such a mechanism depends only on the differ- 
ent thermal expansivities of the R—0 and Ni— O bonds, 
and does not produce softening. The fluctuating tetrago- 
nal/orthorhombic Jahn- Teller distortion of the NiOg oc- 
tahedra necessary to justify the observed jump in the 
elastic modulus is estimated as ^ 3% and the compati- 
bility of such a requirement with other structural studies 
is discussed. 
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distortion. The generic strain s can be decomposed as 
k) = ^7 1^'')' or explicitly £y = ^jT^j with com- 
ponents Ej — {T'^ I e) , the scalar product being defined as 
Ty{T^ : e} — TijEji. Analogously, the uniaxial stress 
along n is decomposed as 



Thanks to the orthonormality of the T'' only the two 
components cr„ and cr^ will couple with A" and A" and 
we write them as 



(TqW, (Ti, = (ToV 



with 



2n33 rt]^]^ 

2 2 
^11 ~ ^22 



(9) 



V2 



The elastic quadrupoles of the two JT states are written 
as 

A" = AT" , A" = AT" 

where A gives the strength of the distortion. 

The application of cr changes their elastic energies by 
d {Eu - Ey) = -vq {Kj - Kj) '^ji = -vqctoX {u - v) 
and their populations. 



A. Appendix A 

We determine the factor that connects the amplitude 
of softening in Eq. ([6|) to the anisotropy of the elastic 
quadrupole tensor of the octahedral distortions AA. This 
can be done by evaluating the relaxation of the reciprocal 
of the Young's modulus, 8E~^ (n), along a generic direc- 
tion h after application of a uniaxial stress u — gq\\\) 
and then performing the angular average of the relax- 
ation strength, iE~^ (n) jE^^ (n), over n, in order to 
obtain an approximation of the polycrystalline average. 
It is convenient to decompose stress and strain into an 
orthonormal basis adapted to the cubic symmetry^ 



1 



rpU 



1 



1 
1 
1 

-10 
0-10 
2 



rjiV 



1 
0-10 




(8) 



where is the dilatation (irreducible representation 
Aig), T" and tetragonal and orthorhombic shears of 
type Eg , and the three shears of type T2g do not enter in 
the calculation, since they are orthogonal to the Eg JT 



and riy = 1 — n„ change as ~ T^T^ {v — u), where 
it is used Act ^ fceT. This change of populations results 
in a change of the anelastic strain 
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4fcBT 



A^ {u - v) [T" - T"] 



The uniaxial component of de^^^ parallel to h is 
(nlde'^^ln) = ^ijUin^de'^f , and using Eq. (0) or 
(n| T" — T" |fi) = (m — w), the relaxation of the Young's 
modulus along n is 



5E- 



CTO 



4fcBT 



A' [u-vf 



A simple formula for the polycrystalline average of the re- 
laxation strength A (fi) = 5E~'^ (n) /E~^ (n) is obtained 
from the Reuss approximation of uniform stress over the 
grainSfSi taking the angular average (A) ~ (^5E~^ (fi)) E 
where the anisotropy of the material is neglected. By 
using Eq. dU and (nf) = i, (nf) = i, {nln]) 
finally obtains 



nhi'j) — one 



(A) 



CVqX^ 

IbkBT' 



■E 
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It is possible to introduce a mixing of the two nor- 
malized symmetry strains, due to hybridisation with 
the 02p orbitals and delocalisation,— by setting A" — 



as 



1 



I) T" + 



and similarly for A" with < / < 

1. In this case, in the above formulas A should be re- 
placed by A (1 — /), but from anelastic measurements it 
is not possible to evaluate A and / independently, and 
we assume f — 0- 



B. Appendix B 

The softening in Eq. ([5]) is valid above the CO or 00 
transition temperature; below that temperature, the res- 
tiffening depends on the precise nature of the transition. 
In the simplest cases of doublet 00,5^ the susceptibility 
below Tc is given by 



1 



X 



fcBTcosh^ (h/kBT) 



where h is the mean field felt by each ion, acting as order 
parameter, and is solution of the self-consistent equation 

h — fceTc tanh (h/k-QT) . 

For CO, consider the Landau expansion in terms of 
powers of the charge fluctuation order parameter, as in 
Ref. [J?! but limiting to the case that only one order 
parameter Q is relevant. Than the free energy is written 



where a = ao(T — Q). Solving the equilibrium condi- 
tions dF/dQ = and dF/dQe = one finds (T) = 
ao//3 (T - Tco) with Tco - + AT, AT = (Coao)- 
The renormalized elastic constant is then 



C = Co-g 



dQ 

de 



where dQ /de can be obtained by derivating with respect 
to e the equation dF/dQ — 0: 

dQ g 



de a + 3^Q2 

When T > 9 it is Q = and we obtain Eq. (O as before, 
here rewritten as 



C^Co 



(T - e - AT) 



(T-e) 

In the CO phase, we obtain 

{T-e- AT) 



C = Co- 



(T-e- |Ar) ■ 

These formulas have been used for the fits in Fig. S) 
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